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In a recent experiment, Segr£ et al. [1] used the particle 
imaging velocimetry (PIV) technique to measure the spatial 
correlation function, ) = (fr(r)^r(r -f 0)> of the veloci- 
ty fluctuation b r in a sedimenting suspension of non-Brownian 
particles over a wide range of particle concentrations and sam- 
ple sizes. They found that the measured C( f) ~ exp(-f/£h 
where the velocity correlation length £ depends on the particle 
radius a and volume fraction do in a non-trivial power-law 
form £ ~ In this paper we propose a new set of 

coarse-grained equations of motion to describe concentration 
and velocity fluctuations in a dilute sedimenting suspension 
of non-Brownian particles. With these equations we find that 
colloidal sedimentation is analogous to high Rayleigh number, 
high Prandtl number turbulent convection [2,3]. Our mod- 
el explains the experimental results by Segit ct al. and also 
provides a coherent framework for the study of sedimentation 
dynamics in different colloidal systems. 

To understand the basic principles governing the colloidal 
sedimentation, we consider a simple case of a dilute sediment- 
ing suspension of hard spheres in a long cylindrical tube of 
radius L. To separate the velocity fluctuation bu from the 
mean settling velocity f, we choose a uniform suspension 
with v = 0 as our reference system. It has been suggested 
[1,4] that velocity fluctuations in a sedimenting suspension 
may arise from fluctuations of the local particle concentra- 
tion. Therefore, we model the colloidal sedimentation with a 
coarse-grained Navier-Stokes equation. The fluid velocity bu 
and pressure bp at a point x satisfy the creeping flow equation 

[5] 

Vbp[x) - ijV 2 $u(x) = f^n(x), (I) 

where // is the viscosity of the fluid and bn [= n(x) — n] 
represents the fluctuation of the particle number density » (x ) 
about its mean ft. In the above, f = (47r/3)n 3 A^g is the 
buoyancy force acting on a particle of radius a, where g is 
the gravitational acceleration and Ap — p v — p* is the density 
difference between the particle ( p p ) and the solvent (/?-). In 
writing Eq. (1) we have assumed that the fluid volume ele- 
ment b\' is a coarse-grained volume, which is large enough 
to contain many particles but is small enough such that the 
particle distribution inside bV is uniform. In this case, we 
have f>n(x) = Apgfa(x) - rio], where <£(x) is the particle 
volume fraction and Oo is its mean value. 

Nondimensionalizing Eq. (1 ) with respect to the length L, 
the time L 1 / D, and the concentration <j> o, we have 

--Vfy(x) -f V^ufx) = Rao(x) z, (2) 

G 


where the unit vector z is directed upward opposite to the 
direction of g, and the dynamic pressure bp has included a 
term, -Apg&oz, to absorb contributions from the constant 
forcing term — A/?g<3o* In Eq. (2) the Rayleigh number Ra is 
defined as 

Ra = Apg<faL 3 f(yD), (3) 

where D is an effective diffusion constant of the particles. 
The Schmidt number a is given by <t = vj D with v being 
the kinematic viscosity of the fluid. For a dilute suspension 
of small colloidal particles, D is approximately equal to the 
particle self diffusion constant = IbT /( toga), where 
k B T is the thermal energy. For large non-Brownian particles, 
however, the effect of thermal agitations is negligible and their 
diffusion-like motion is produced by the hydrodynamic inter- 
actions between the particles [6]. Nicolai etal. have shown [7] 
that the hydrodynamic diffusivity has the form Dh — 5«f : a» 
where f r o = 2a 2 Apg/{9g) is the Stokes velocity. 

With the hydrodynamic diffusivity Dh, Eq. (3) becomes 

Ra =0.9o 0 (-) 3 - (4) 

a 

It should be mentioned that while it is cancelled out in Ra, 
A py is needed so that D fl can be used to describe the hy- 
drodynamic diffusion of the settling particles at small length 
scales. Equation (2) together with the continuity equation for 
an incompressible fluid 

V-*u = 0 (5) 

and the adveclive mass diffusion equation 

d t 6 + [Hi- V)<? = V 2 <? (6) 

complete the description of concentration and velocity fluctu- 
ations in colloidal sedimentation. 

It is evident that Eqs. (2)-(6) are the same as those for 
buoyancy -driven convection [3]. Velocity and concentration 
fluctuations in colloidal sedimentation arc therefore analogous 
to those in buoyancy-driven convection, and they are com- 
pletely controlled by the two dimensionless parameters Ra 
and g , once the boundary conditions arc specified. We now 
estimate typical values of Ra and o in colloidal sedimenta- 
tion. In the experiment by Segr£ et al. [1], the particle’s radius 
a ^ 8 pm, Stokes velocity Uq ~ 6.5 pm/s, volume fraction 
(f)Q cz 0.05, and the characteristic sample size L ~ 1 cm. With 
these experimental values, we find Dh . — 2.6 x 10 6 cm 2 /*, 
Ra ~ 8.8 x 10 7 , and g ~ 3800. The Schmidt number a 
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is equivalent to the Prandtl number in thermal convection. 
Colloidal sedimentation is, therefore, associated with high 
Rayleigh number, high Prandtl number turbulent convection. 

To understand the sedimentation dynamics, it is helpful 
to distinguish two characteristic length scales in convection: 
the viscous dissipation length b v and the diffusive dissipation 
length b d . The values of b v and b d are determined, respective- 
ly, by the transition Reynolds number Re c = bubv/ v and the 
transition Peclel number Pe c = bubal Dh • Here bu is the rms 
value of the velocity fluctuation bu averaged over a volume 
of bl (or 6d 3 ). It is the ratios of these lengths to each other 
and to the sample size I that determine the flow state of the 
system [2]. For high -Ra, high -rr turbulent convection, one 
anticipates that the flow consists of three different regions: (i) 
a < ( < b<u (ii) b d < C < b l t and (iii) b d < C < L. In 
Region (i), molecular viscosity and hydrodynamic diffusivity 
determine the momentum and mass transport processes, re- 
spectively, and hence the particle distribution remains uniform 
without any large fluctuations. In Region (ii) turbulent (or 
eddy) diffusivity and molecular viscosity are dominant, and 
thus large fluctuations in particle concentration are expected 
but the velocity field remains relatively smooth. Finally, in 
Region (iii) turbulent diffusivity and viscosity both dominate 
over the corresponding hydrodynamic and molecular process- 
es. In this case, one expects to see large fluctuations both in 
particle concentration and in velocity at different length scales. 

We first discuss the length b d , above which velocities 
become large and concentration fluctuations are transported 
by convection. This occurs when the local Peclet number 
Pt ~ bu(/Dh becomes larger than Pe c . Recent thermal 
convection experiments have shown [8] that while turbulent 
mixing creates on average an isothermal fluid in the turbulent 
bulk region, large temperature fluctuations still remain in the 
region and the characteristic length scale associated with these 
fluctuations is of the order of b d . Therefore, the velocity cor- 
relation length £ is determined by b d in Region (ii). According 
to Kraichnan’s theory [21, 

bj- (2x 2 Pc l c )' 13 LRa-'* 3 , (7) 

where the power law amplitude is expressed in terms of the 
numerical value of Pc L . Priestley [9] first gave a direct argu- 
ment for the Ra~ l/2 scaling. He argued that when Ra is large 


enough, b d should be a new length scale independent of the 
sample size L. With Eqs. (7) and (4), we immediately have 
£ ~ b d ~ LRa ~ l/ 3 ~ a<f>~ U \ The mapping of colloidal 
sedimentation to turbulent convection, therefore, explains the 
experimental finding that £ ~ 1 1 a.<6 0 Ii provides a 

physical interpretation for the existence of a velocity cut-off 
length, which prevents hydrodynamic dispersion coefficients 
from being divergent. 

We now discuss the velocity variance b u in colloidal sed- 
imentation. According to Kraichnan’s theory [2], 

~ Pt c Dh PtcDh /o\ 

bu — — — 7773 ■ 1°) 

(27T ipe 2 c) 1 LRa-W 

Eq. (8) states that at the transition Peclel number Pc c , the 
mass flux due to hydrodynamic diffusion, Dhb$ j b d> is ap- 
proximately equal to that by convection, bubtp. Because 
b d ~ O0Q 1/3 and D h ~ «Fo, we find from Eq. (8) that 
bv ~ bu ~ Dh/bd ^ DoC>l n , which is independent of the 
sample size L. This result agrees well with the experimental 
finding that b v — 21 [1]. 
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